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{A^) Condensate, Bianchi Identities and Chromomagnetic 
Fields Degeneracy in SU(2) YM Theory 

F.V. Gubare^O and S.M. Morozo\0 
Institute of Theoretical and Experimental Physies, 

B.Cheremushkinskaya 25, Moscow, 117259, Russia 

We consider the non-Abelian Bianchi identities in SU(2) pure Yang-Mills theory in D=3,4 focusing on 
the possibility of their violation and the significance of the chromomagnetic fields degeneracy points. 
We show that the recently proposed non-Abelian Stokes theorem allows to formulate the Bianchi 
identities in terms of the physical fluxes and their relative color orientations. Then the violation 
of Bianchi identities becomes a well defined concept ultimately related to the degeneracy points. 
The locality and gauge invariance of our approach allows to study the problem numerically. We 
present evidences that in D=4 the suppression of the Bianchi identities violation is likely to destroy 
confinement while the removal of the degeneracy points drives the theory to the topologically non¬ 
trivial sector. However, confronting the results obtained in three and four dimensions we argue that 
it is the mass dimension two condensate which probably explains our findings. 
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I. INTRODUCTION 

Gauge theories are usually formulated in terms of the 
gauge potentials taking values in the Lie algebra 
of the corresponding gauge group. Provided that the 
gauge coupling is small this description is indeed ad¬ 
equate and provides local functionally independent co¬ 
ordinates on the configuration space. However, in the 
strongly coupled gauge theories the potentials themselves 
obtain a separate physical meaning. Here we mean the 
non-perturbative dimension 2 condensate intro¬ 

duced in HQ, which received a particular attention in 
recent years (see, e.g. Ref. for review and further 
references). 

The original motivation of this work was the analysis of 
various possible contributions to the condensate. 

Note that the central point of Ref. [y was in fact the 
consideration of the Abelian Bianchi identities and their 
ultimate relation to As far as the Abelian the¬ 

ory is concerned the non-triviality of {A^^A condensate 
is essentially equivalent to the Bianchi identities viola¬ 
tion. Therefore in the non-Abelian case it seems natural 
to start from the corresponding Bianchi identities and 
investigate their role in the condensate forma¬ 

tion. However, the literature on the subject turns out 
to be scarce. In particular, as is well known from the 
Abelian models the rigorous treatment of the Bianchi 
identities requires the non-perturbative (say, lattice) reg¬ 
ularization. But we were unable to find papers devoted 
to this problem in the non-Abelian case. 

On the other hand, the investigation of the non- 
Abelian Bianchi identities is important on its own right. 
Without mentioning all the aspects of the problem, let us 
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note that the (^^m) condensate is certainly connected 
with the non-Abelian Bianchi identities. Moreover, it 
was emphasized in Refs. Q, Q Q that the Bianchi iden¬ 
tities and the possibility of their violation are ultimately 
related to the confinement problem. Then the logic sug¬ 
gests to consider whether the (^^m) condensate is rel¬ 
evant for confinement as well, the question which was 
discussed in Qi (see also 0)- Therefore we see that all 
these problems are in fact indispensable from each other 
and cannot be considered separately. We decided to fo¬ 
cus on the Bianchi identities in this paper; the connection 
with the quantities like is discussed in the due 

course. Throughout the paper we work with Euclidean 
three and four dimensional SU(2) gluodynamics keeping 
in mind the lattice regularization of the theory, although 
we nowhere relay exclusively on the lattice. The paper is 
reasonably self-contained, the results which we’re using 
are briefly reviewed. Note that the similar in spirit but in 
no way identical treatment could be found in Refs. HQ 

The primary tool of our analysis is the non-Abelian 
Stokes theorem [13 derived recently by one of us. The 
advantage is that it allows to work directly in terms of 
the gauge invariant quantities like magnitudes of the ele¬ 
mentary fluxes and their relative orientations. As might 
be expected the non-Abelian Bianchi identities could be 
reduced to the application of the above theorem to the 
infinitesimal closed surfaces. However, in this case the 
non-Abelian Stokes theorem not necessary gives zero, the 
answer, in fact, is proportional to the integer number. 
Since every step in the derivation is gauge invariant this 
integer is gauge invariant as well and in the continuum 
language corresponds to the non-Abelian Bianchi identi¬ 
ties violation. 

The non-Abelian nature of the theory manifests itself 
in the complicated geometry underlying the Bianchi iden¬ 
tities. We consider all these questions in detail and show 
that the careful but purely geometrical treatment leads 
to the consideration of the special degenerate points in 
the configuration space at which a particular determi- 
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nants constructed from chromoelectric and chromomag- 
netic fields vanish. Finally we show that the investiga¬ 
tion of the non-Abelian Bianchi identities is indispensable 
from the study of these degenerate points. Therefore the 
framework outlined above naturally extends to include 
the degeneracy points, the relevance of which both for 
confinement and chiral symmetry breaking was discussed 
in Refs. Him. 

The locality and gauge invariance of our construction 
allow us to study the problem numerically. We inves¬ 
tigate the effects due to the Bianchi identities violation 
and the degenerate points in the numerical simulations. 
As might be a priori expected the suppression of the 
degenerate points always lead to the violation of the re¬ 
flection positivity. Moreover, in D=4 one could easily 
pin-point the origin of the reflection positivity violation: 
it is caused by rapidly rising global topological charge. 
Thus in D=4 the suppression of the degenerate points 
shifts the vacuum to the non-trivial topological sector. 

As far as the Bianchi identities are concerned the re¬ 
sults depend crucially on the space-time dimensionality. 
In D=3 the suppression of the Bianchi identities violation 
does not change the theory in any notable way. How¬ 
ever, in D=4 the effect is different: it seems that the 
suppression of the Bianchi identities violation is likely 
to destroy confinement while other measured character¬ 
istics of the theory remain qualitatively unchanged. At 
least this is so for the lattices and coupling constants we 
have considered. Note that the problem still requires a 
careful numerical investigation, in particular, we had not 
studied yet the volume dependence of our results. The 
corresponding analysis will be published elsewhere. 

Finally we argue that it would be misleading to in¬ 
terpret our results as the statement that confinement is 
caused by the Bianchi identities violation. Confronting 
the results obtained in three and four dimensions we show 
that it is the condensate which is probably rele¬ 

vant for confinement. Although the argumentation is not 
rigorous it seems to be the only one which matches our 
findings. 


II. FORMULATION OF THE PROBLEM 


The primary object of our investigation is the Bianchi 
identities for SU(2) gauge fields in four space-time di¬ 
mensions. Thus we will analyze the equations 

( 1 ) 

having in mind eventually Euclidean lattice regulariza¬ 
tion of SU(2) pure Yang-Mills theory. Here E®,, is the 
conventional continuum field-strength tensor 

F;^ = dpAl-d,Al + e<^^‘^AlAl, ( 2 ) 


Greek and Latin indexes run through 0, ..,3 and 1, ...,3 
respectively. Our treatment also applies in three dimen¬ 
sions where Bianchi identities are as follows 


d, Bf + = 0 , 

= [D = 3]. 


(3) 


However, it turns out that the three-dimensional case is 
physically quite different from D = 4 and we’ll comment 
on that in the due course. 

In this section we give qualitative continuum argu¬ 
ments which show that at least at some points in the 
configuration space the Bianchi identities should 

be considered with care. 


A. Chromomagnetic fields degeneracy 

It has been known for a long time that in non-Abelian 
gauge theories two or more gauge inequivalent potentials 
could produce the same field strength M This phe¬ 
nomenon, known as Wu-Yang ambigui^ had received 
great attention in the past (see, e^. [1^ [Tj, [I^) and 
it was noted long ago [3 13 Fml that in D=4 the 
Bianchi identities constitute an algebraic obstruction for 
the ambiguity to exist. Namely, for given chromoelectric 
Ef = Fgj and chromomagnetic Bf = 1/2 Sijk F^f. fields 
Eq.® is a linear algebraic system of 12 equations for 
12 unknown A“. Therefore away from the set of points 

where the matrix F/i^^ degenerates 

detT = 0 (4) 

Bianchi identities allow to express the gauge potentials as 
local single-valued function of E/ and B/. On the other 
hand, there is no physical principle or symmetry which 
could keep the sign of det T fixed. Indeed, in the weak 
coupling perturbation theory the sign of det T changes 
wildly and therefore the degeneracy of chromomagnetic 
fields, Eq. ®, is, in a sense, generic. Note that the sit¬ 
uation is quite different in D=3 since Eq. m formally 
constitutes 3 equations for 9 unknown variables. There¬ 
fore in three dimensions the Bianchi identities do not con¬ 
strain the gauge potentials at all and the Wu-Yan g am bi- 
guity problem is much more severe (see, e.g. Refs. flHl^ 
for discussion). Unfortunately, we are not aware of any 
conclusive considerations of the degenerate points 0 in 
the literature. It is true that Eq. 0 by itself is known for 
a long time but most of the analysis performed 

so far considered it in the context of dual formulation of 
gluodynamics El mil in 123 from which the informa¬ 
tion about original Yang-Mills fields is hard to extract. 
Ref. El seems to be the only exception where it was ar¬ 
gued that physical wave functionals should vanish at the 
points of degeneracy. We will see below that equations 
similar to 0 arise naturally in the construction of the 
Bianchi identities. Moreover, the points of degeneracy 
seem to be relevant for gauge fields dynamics. 
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What we have said so far is in accordance with gen¬ 
eral expectation that in the non-Abelian gauge theories 
there is no unique way to express in terms of the 
corresponding field strength (apart from the usual gauge 
ambiguity, of course). At this point one could give an ex¬ 
ample of special gauges (complete axial, coordinate, con¬ 
tour gauges, see for review), in which the gauge 

potentials are always explicit single-valued functions of 
the field strength. Is there any contradiction? Although 
this question is not directly related to our work, we note 
that all the gauges mentioned above are consistent only if 
Bianchi identities o, ® are satisfied identically [sj • In 
particular, in the Abelian case one notices that the 
presence of elementary magnetic charges forces the po¬ 
tentials in contour gauge to depend upon the arbitrary 
contour prescription. Of course, this is a manifestation of 
famous Wu-Yang ambiguity which in this case certainly 
arises because point-like monopoles violate the Bianchi 
identities. We conclude therefore that the possibility of 
Bianchi identities violation should not be excluded a pri¬ 
ori. Moreover, the very existence of Wu-Yang ambiguous 
potentials hints on the violation of CJ, ®. 


B. Bianchi identities violation 

The possibility that the r.h.s. of Eqs. o, © might be 
non-zero was considered long ago (see, e.g. j^), but as 
far as we know this approach had never been actively 
developed. This is mostly because the study of Bianchi 
identities violation requires a particular regularization, 
which should correctly respect the global structure of the 
gauge group. It turns out that for our purposes the lat¬ 
tice formulation is distinguished (see Refs, H^l for 
discussion). Therefore consider the basic SU(2) gauge 
theory observable, which is also the fundamental object 
on the lattice, the Wilson loop in spin 1/2 representation 


flux is always bounded 0 < $(C') < tt due to periodic¬ 
ity (compactness) of the gauge action. Moreover, there 
exist no physically meaningful experiment which could 
distinguish the fluxes <&(C') and $((7) -I- 27r and this ob¬ 
servation applies equally well to the infinitesimal con¬ 
tours which constitute the lattice definition of the field 
strength. On the other hand, there is no trace whatso¬ 
ever of the gauge action compactness in the continuum 
expression ©■ In this respect the SU(2) gluodynam- 
ics is similar to the compact U(l) gauge model [SJI (see 
Ref. for review). In fact, some consequences of the 
compactness of the non-Abelian gauge theories were al¬ 
ready discussed in the past [s^. Note however that we 
are not saying that singular fluxes are important in the 
continuum limit of lattice formulation. After all this is 
a dynamical question which cannot be studied with sim¬ 
ple arguments above. Rather we would like to point out 
that the very definition of on the lattice is a priori 
different from the continuum one © and therefore the 
validity of O, © in the lattice context should be con¬ 
sidered anew. We stress that our arguments are purely 
kinematical and follow directly from the gauge invariance 
along. Whether or not the violation of Bianchi identities 
is physically relevant is a dynamical issue which we in¬ 
vestigate (at least partially) later on. 

To conclude we note that nowadays there exist both 
theoretical arguments iii and the experimental lat¬ 
tice data which favor the non-vanishing r.h.s. of 
Eqs. o> © in the continuum limit of lattice gauge mod¬ 
els. Although the approaches of these papers are quite 
different, the conclusion is essentially the same: the non- 
Abelian Bianchi identities are indeed violated in the scal¬ 
ing (continuum) limit and this fact is related to the prob¬ 
lem of confinement. 


III. LATTICE BIANCHI IDENTITIES 


W(C, Xo) = Pexp icr“ A“ dxi ^, 

W{C) = ^TiW{C,Xo). 


(5) 


Here tT“ are the Pauli matrices, C is some closed con¬ 
tour with marked point xq € C from which the path or¬ 
dered integral starts and P-ordering is defined from left 
to the right. Note the unusual normalization of SU(2) 
generators which we take for future convenience. By def¬ 
inition the operator W{C,Xq) measures the non-Abelian 
flux <i)(C, Xo) penetrating the contour 


W{C, Xo) = -J-YCxo) ^ ^(c') = cos $(C), 

( 6 ) 

$(C) = v'$“(C',xo)$“(C,xo), 

where the flux $(1^) is gauge invariant and does 
not depend on xq- Eq. © will be thoroughly analyzed 
later, but now we note that the physically observable 


A. Preliminaries 

In this section we briefly summarize what has been known 
so far about the non-Abelian Bianchi identities on the 
lattice and comment on the strategy we employ in this 
paper. Surprisingly enough the literature on the subject 
seems to be very scarce (contrary to the Abelian case 
which we do not consider however) and the most rele¬ 
vant for our discussion references are |3^ (see 
also [ 3 ^). Historically, the Bianchi identities explicitly 
appeared first in the context of plaquette (field-strength) 
formulation of lattice QCD |3^ Idfl . In particular, it 
was noted that the strong coupling expansion can be 
obtained as an expansion towards restoring the lattice 
Bianchi identities. 

It turns out that the formulation of Ref. is the 
most appropriate for our purposes. Essentially it con¬ 
sists in the observation that any lattice gauge field con¬ 
figuration could be interpreted as a homomorphism from 
the lattice edge path group into the gauge group (see 



Ref. EH for definitions). It follows form the definition of 
homomorphic mapping that 
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U{C,y o C-1) = U{C,y)U-\C,y) = 1, (7) 

where C^y is arbitrary path connecting the points x and y 
and the composite path CxyoC~y is usually referred to as 
null-homotopic. In fact, Eq. o looks rather obvious for 
everyone familiar with lattice formulation. However, the 
assertion of Ref. is that Eq. 0 constitutes the most 
general form of lattice Bianchi identities and indeed just 
that: an identity. Note that Eq. o looks quite different 
from what is expected in the continuum. To establish the 
relation between o and 0-0 consider the path C^x 
shown on Figure ^ It follows trivially that the equality 
U{Cxx ° Cxx) = 1 is equivalent to 

U{Ri)U{R2)U{R^)U{Ri)U{R^)U{R&) = 1, (8) 

which in the naive continuum limit reduces to the con¬ 
ventional Bianchi identities 1^, Q. Moreover, Eq. © is 
the particular case of the so called operator non-Abelian 
Stokes theorem (see, e.g. for review) which 

allows to represent (rather formally though) the path or¬ 
dered exponent as the surface ordered integral 

Pexp i J dx^ = Ps exp ^ J , (9) 

C=SSc Sc 

where ^ is non-local covariantly transformed field- 
strength the concrete form of which is not important for 
what follows. The surface Sc is arbitrary and consis¬ 
tency requires the representation 0 to be independent 
on Sc as long as SSc = C. In particular, the r.h.s. 
of Eq. @ being applied to closed surface Sq, SSq = 0, 
should always give the identity 

Psexp ^ f = 1. (10) 

^ J Sq ,(5So—0 

In fact, Eq. ® is the special case of 113 in which Sq is the 
boundary of elementary lattice cube. Therefore, it seems 
to be legitimate to formulate the non-Abelian Bianchi 
identities as the requirement of surface independence of 
the non-Abelian Stokes theorem. 

Eqs. 0-1101) are the starting point of our considera¬ 
tions below. However, before going into details let us 
comment a bit on our strategy. We note first that the 
identity on the r.h.s. of Eqs. ®, cnj) could in general be 
written as 

1 = = 1, q € Z. (11) 

The color direction ft is gauge variant and will not con¬ 
cern us here. Suppose that we are able to give an un¬ 
ambiguous gauge invariant meaning to the integer q and 
that it is non-zero for some Sq in given gauge background. 
Then this would certainly mean that there is a point 
somewhere inside Sq at which the continuum Bianchi 
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FIG. 1: Graphical representation of lattice Bianchi identities. 

identities are violated. Here the argumentation is essen¬ 
tially the same as in well known Abelian case. So the 
problem is to make sense of q which should be well de¬ 
fined and gauge invariant. From now on we refer to the 
integer q as the “magnetic charge” whatever it is. In par¬ 
ticular, neither charge conservation nor any other usual 
properties of the magnetic charge are assumed. Secondly, 
Eqs. 0,03) are not quite suitable to analyze the Bianchi 
identities. This is precisely because neither 0 nor ca 
make, in fact, no direct reference to the non-Abelian field 
strength. And this is in sharp contrast with the Abelian 
theory in which the Bianchi identities even on the lattice 
explicitly refer to physical fluxes. It turns out that the 
solution of the second problem simultaneously solves the 
first, namely, the non-Abelian Stokes theorem being ex¬ 
pressed in terms of the physical field strength provides 
the definition of q which we are looking for. 

B. Chromomagnetic Fields on the Lattice 

The distinguished feature of the lattice regularization is 
that the gauge theory is formulated in terms of the Wil¬ 
son loops along and strictly speaking the lattice does not 
need to introduce the notion of the field strength. Chro¬ 
momagnetic fields appear only in the limit of vanishing 
lattice spacing, overwise one should rather think in terms 
of the non-Abelian fluxes which are defined by Eqs. ©, 
0. Therefore consider the Wilson loop 

W{CA) = Pexp A(T)dr = , 

A{T)=a‘^A-{x)x>^{T), n^{C,t) = l, (12) 

W{C)= iTrW(C',t) =cos d>(C'), 

defined for some closed contour C = {a:(t),0 < t < 
T, a:(0) = x{T)} (our presentation is similar but not 
identical to that of Ref. [Tnji . see also Ref. 0). We 
assume that W{C) ^ ±I and then it is convenient to 
parametrize the Wilson loop in terms of the flux magni¬ 
tude $((7) G (0; tt) and the instantaneous flux direction 
in color space n{C, t) which explicitly depends on t. It is 
clear that $((7) is gauge invariant while n{C,t) rotates 
as three-dimensional vector under the gauge transforma¬ 
tions at point x{t). Consider now another contour C 
which touches (or intersects) (7 at point x{to) = 
Evidently, while both n{C,to) and n(C",tg) are gauge 
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variant their relative orientation (angle in between) is 
gauge independent. Moreover, the construction could be 
iterated: for any number of contours intersecting at one 
point the relative orientation of instantaneous fluxes at 
that point is gauge invariant. It is amusing to note that 
the relative orientation of elementary fluxes received al¬ 
most no attention in the past. While the magnitude of 
various fluxes had been discussed and measured in var¬ 
ious circumstances (see, e.g. Ref. and references 
therein), it seems that only Refs. [iM 14^ studied their 
relative orientations. 

Consider next the behavior of the flux parametrized by 
Eq. (I12II under the change of contour orientation. Phys¬ 
ically one expects that the total flux should change sign 
when contour is followed in the opposite direction 


it is natural to describe the instantaneous flux direction 
by means of fictitious (iso-)spin 1/2 particle living on 
the contour. The spinor wave function is given by two- 
component normalized complex quantity 

{z\ = [z^-Z2], {Z\z) = \z,\’^ + \Z2\^ = 1, (16) 

which is bra-vector in accordance with our left to the 
right P-ordering convention. The defining equation for 
the Wilson loop becomes the Schrodinger equation for 
spinor 


(z(t) \{idt + A) = 0, 

( zit) I = (^(0) I • Pexp i J* A(t) dr . 


(17) 


$“((7-1, t) = $(C'-i)n“(C'-\t) = -$“(C',t). (13) 

The parametrization m respects the intuition and in¬ 
deed the flux direction changes sign while the flux mag¬ 
nitude is orientation independent 

n(C-\t) = -n(C,t), ^iC-^) = ^{C). (14) 

Here we come to the important point concerning the 
determination of physical field strength from the in¬ 
finitesimal fluxes. Suppose that we measure twice the 
elementary flux, first with an oriented area element Scr^‘^ 
and then with reversed orientation 6a''^ = Evi¬ 

dently, the corresponding Wilson loops are conjugated to 
each other 


Therefore the Wilson loop eg is the quantum mechani¬ 
cal evolution operator for spin degrees of freedom. As is 
usual in quantum mechanics the state vectors could be 
arbitrary rephased 

{z{t) \ ^ e*^(‘)(z(t)|. (18) 

The particular choices Imzi = 0, lmz 2 = 0 lead to 
well known families of (anti)holomorphic spin coherent 
states [53 (see, e.g. [HJ for review). Following the quan¬ 
tum mechanical analogy 113,13 could argue that the 
eigenstate of the evolution operator W (C, 0) 

(z(0)|W(C',0) = e*^(^)(z(0)|, (19) 


lT(fc'^") = W\Sa''>^). (15) 

On the other hand, the expansion in powers of lattice 
spacing a reads 

W{6a>^‘') = 1 -b + 0(a'‘), 

= 1 -b id+ 0{a*) = W{Sa^''). 

and disagrees with lisj. This simple exercise which ap¬ 
plies equally in the Abelian case shows that the lattice 
area element dx^dx'^ is in fact unoriented dx^dx'' = 
dx'^dx^ contrary to the usual continuum relation = 
dx^ A dx'^ = —dx'^ A dx^. Therefore in order to define 
the field strength on the lattice a canonical orientation of 
all elementary squares (plaquettes) should be fixed first. 
Overwise the field strength will suffer from sign ambigu¬ 
ity on different plaquettes. In fact, the canonical ordering 
is well known in lattice community and the conventional 
agreement is to consider with ^ < v only. However, 
the orientation conventions are crucial for the interpre¬ 
tation of lattice equations below in the continuum terms. 
From now on we always assume that the infinitesimal 
fluxes are constructed with canonically oriented plaque¬ 
ttes. 

It is convenient to generalize the representation 03 in 
order to gain a simple physical interpretation. Namely, 


is of special importance and is usually referred to as cyclic 
state. In particular, the state (z(0) | being the eigenstate 
of 1E(C', 0) at t = 0 remains the eigenstate of W{C,t) 
during the evolution It follows immediately that 

the cyclic state (HU is best suited to describe the instan¬ 
taneous flux direction. Indeed, it is a matter of one-line 
calculation to show that n°'{C,t) = {zit)\a°‘ \z(t)). 

In other words, the flux direction n(C, t) and the ratio 
Z 2 it)/zi{t) of cyclic state components are related to each 
other by standard stereographic projection. In particu¬ 
lar, the flux magnitude is given by 

<i>(C) = arg[(z(t)|W(C,t)|z(t))] (20) 

and is t-independent. Moreover, if contour C is subdi¬ 
vided into N segments then 


$(C) 


N-1 

arg (z(tfc) |Pexp 


ife+l 

iJM t) dr I z{tk+i)) , 

tk 


( 21 ) 


where the identification to = tjv is assumed. As far as the 
relative orientation of fluxes is concerned it is tempting 
to consider the quantities like arg(z|C)- However, it is 
not invariant under m because ( z | and ( C | could be 
rephased independently. Nevertheless, the equations we 
will get do indeed include the products like ( 2 ; | C) yot 
respecting the U(l) invariance (ITHll . 
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It remains only to consider the multivaluedness of the 
cyclic state defining equation m- Indeed, there exist 
two solutions of Eq. (P|l while we discussed only one of 
them. The second eigenstate is obtainable from the first 
one by substitution 

^2 ^ zl zi -z^. (22) 

It is clear that Eq. (Unj corresponds to the “spin-up” wave 
function for which the spin is aligned with the magnetic 
field, while the second eigenstate (1^ is the “spin-down” 
state which has spin anti-aligned. Our original goal was 
to describe the direction of instantaneous flux and there¬ 
fore the anti-aligned state should be discarded since it 
corresponds to the inverted flux direction. Note also 
that the flux magnitude $(0) is positive by definition 
but with anti-aligned state we get d)(C') < 0. We con¬ 
clude therefore that for given contour orientation there 
is no ambiguity in Eq. (PI and the appropriate family 
of cyclic states {z{t)\ is uniquely defined. The second 
“spin-down” eigenstate describes the flux direction for 
inverted contour orientation and therefore Eq. ipil cor¬ 
responds to the time reversal operation for spinors in 
quantum mechanics. 

The above considerations apply immediately on the 
lattice. The only difference with the continuum is that 
the gauge potentials are unknown, we have only the par¬ 
allel transporters along the elementary links. But this is 
actually enough: the Wilson loop is constructed by direct 
matrix multiplication and then Eq. m applies literally. 
The instantaneous flux direction is determined via m 
or 03 at lattice sites passed by Wilson loop. The flux 
magnitude is given by Eqs. (P|. (ET| 

To summarize, every Wilson loop W (C, t) is char¬ 

acterized by the magnitude of the flux $((7) and the in¬ 
stantaneous flux direction n[C,t), = 1 which varies 

along the contour and is reversed on changing contour 
orientation. The quantum mechanical language is ade¬ 
quate to describe both $((7) and n{C,t): there is a fic¬ 
titious spin 1/2 particle living on (7, the polarization of 
which gives exactly n((7, t); the wave function of the par¬ 
ticle is defined for given gauge background uniquely up 
to the phase and change of contour orientation is equiv¬ 
alent to time reversal operator applied to the spinor; the 
particle evolution along C is cyclic, initial and final states 
differ only by phase and this phase is the magnitude of 
the flux penetrating C. On the lattice the difference is 
that the flux direction (wave function of spinning parti¬ 
cle) is known only at lattice sites x G C. Moreover, the 
orientation of all elementary plaquettes is fixed to be the 
canonical one. 


C. Non-Abelian Stokes Theorem 

The last ingredient which we need to complete the pro¬ 
gram outlined in sec. lIII Al ls the non-Abelian Stokes the¬ 
orem derived recently by one of us M- Although the 
results of Ref. Hi are applicable almost literally, let us 
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FIG. 2: Segment of the Wilson loop W(C) in the original 
and ribbon-like representation. The operator in between solid 
blobs is Uk,k+i, Eq. 11281 . 


review them in order to introduce the notations and com¬ 
ment on the differences with present work. 

Therefore consider the Wilson loop W{C), segment of 
which is shown by straight horizontal line on Figure |21 
and the surface Sc bounded by (7, which is to the top 
of contour on the same figure. According to what had 
been said above we assign to every plaquette p G Sc 
and Wilson loop itself the corresponding flux magni¬ 
tudes $(p), d)((7) and the instantaneous flux directions 
{zk{p) I, (Cfe(C) I correspondingly (plaquette vertices are 
followed according to the orientation induced by C while 
the states (zu (p) \ are constructed in accordance with the 
canonical orientation). It is convenient to use the graphi¬ 
cal ribbon-like representation in which all plaquettes and 
Wilson loop contour are slided apart (Figure [2|. Let us 
denote 


Uk,k+i = Pexp i f A{t) dr (23) 

and consider the matrix element 

{CkiC)\Uk,k+i\Ck+i{C)) = const (24) 

where const is some real positive number which is irrel¬ 
evant. According to 


$(C) 


Y.cl,k,k+i{C) 

_ k 


mod 27t . 


(25) 


The important observation of Ref. M is that the ma¬ 
trix element (1^ might be calculated in (Zk (p) \ basis 
provided that the relative orientation of plaquette and 
Wilson loop fluxes is taken into account 


{(:k{C)\Uk,k+i\(:k+i{C)) = const ■ {a{C)\zk{p)) X 
X {zk{p)\Uk,k+i\zk+i{p)){zk+i{p)\Ck+iiC)). (26) 

The equality (PI was shown in Ref. M in the matrix 
form. Here we note that Eq. (uni follows from its invari¬ 
ance under m and the unitarity of the evolution oper¬ 
ator (P|. In fact, the relations similar to (Pl are well 
known in quantum mechanics 1521 (see, e.g. |53L l54L 
for details). In particular, Refs. showed the im¬ 

portance and ph ysic al significance of the geodesic inter¬ 
polation used in 
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FIG. 3: fix (left) and '^x (right) terms in Eq. 11271 . Arrows 
correspond to the orientation induced by C. 

Applying Eqs. ((23), ® repeatedly for every link 
of Sc one gets the non-Abelian Stokes theorem 

<i>(C') = ^ I{p) $(p) + ^ ^ y, + 2Trk{Sc) , 

P^Sc x^Sc x^C 

k{Sc)GZ, (27) 

where $(p) is the plaquette flux, 1/2 Tr VF (p) = cos $(p), 
and the factors I{p) = ±1 are analogous to the usual inci¬ 
dence numbers in the differential geometry Ell: I{p) = 1 
if vertices of the plaquette p are followed in the canonical 
order and I{p) = — 1 over wise. The remaining terms are 
illustrated on Figure|31 In particular, 

fix = arg[( zi\z 2 ){z 2 \z 3 ){z 3 \za){z 4 ,\zi)] mod 27r (28) 

is the oriented area of spherical quadrilateral polygon [t^ 
(solid angle) in between the flux directions on the pla- 
quettes pi, ...,p 4 . It is known in quantum mechanics as 
Bargmann invariant [5(i| for the particle’s wave functions 
(see, e.g. El El for review). Physically fl^ accounts 
for the difference of flux orientations on the plaquettes 
sharing the same point x. The third term 

lx = arg[(C|zi)(zi |z2)(z2|C)]mod27r (29) 

equals to the oriented area of spherical triangle con¬ 
structed from the Wilson loop flux direction at x and 
the flux orientations of two plaquettes pi,P 2 G Sc touch¬ 
ing C and sharing the point x. Eq. (PI is again the 
Bargmann invariant for the wave functions of three par¬ 
ticles living on C, pi and p 2 - 

Note that we have omitted the mod 27r operation on 
the r.h.s. of Eq. PI and wrote instead the additional 
2Trk{Sc) term, such that <i>(C') G (0; tt). It is clear that 
k{Sc) is not vanishing in general and is analogous to the 
Dirac string contribution in the Abelian Stokes theorem 
applied for compact U(l) gauge fields |3E3 (see EEIH 
for review and further references). This is in accordance 
with the discussion in sec. EH where we noted that the 
SU(2) gauge model is intrinsically compact and is simi¬ 
lar to compact photodynamics in this respect. However, 
in the non-Abelian case the non-zero k{Sc) could come 
from either of three terms in Eq. PI. In particular, the 
Dirac string contribution k{Sc) ^ 0 does not necessary 
corresponds to the singular elementary non-Abelian flux 
(singular field strength). It could equally come from fl^, 

')x terms which are genuine non-Abelian contributions. 

Note that Eq. (12711 is not only invariant under SU(2) 
gauge transformations, it also remains intact with re¬ 
spect to local (gauge) rephasing (Hl^ (this U(I) gauge 


symmetry is crucial for the dual representation consid¬ 
ered in Ref. E3)- We are in haste to add however, that 
this does not concern the 2'Kk{Sc) term. As might be 
expected the Dirac string contribution is not invariant 
with respect to either of the symmetries. Eq. PI could 
be illustrated nicely in the particular case of pure Abelian 
gauge background. In the Abelian limit all fluxes become 
aligned, but their directions could be opposite. For anti¬ 
aligned flux directions the Bargmann invariants (PI,(PI 
become strictly speaking undefined. For instance, the 
area of the spherical triangle (1^ is undefined when two 
of its vertices are at the north pole of the two-dimensional 
sphere while the third one is at the south pole. However, 
we could avoid this degenerate case by changing simul¬ 
taneously the sign of both n{p,t) and $(p) which does 
not affect the parametrization P|. The flux magnitude 
becomes not positively definite and the incidence coef¬ 
ficients could be absorbed into the definition of $(p). 
Then the second and third terms, which account for the 
flux rotation in color space, vanish and Eq. m becomes 
identical to the usual Abelian Stokes theorem. 

To summarize, the flux <&(C') could be represented al¬ 
most entirely in terms of local physically observable con¬ 
tributions coming from the arbitrary surface Sc bounded 
by C. The point of crucial importance is that all these 
terms are “almost total differentials”: without mod 27r 
operation both the plaquette flux (PI and the Bargmann 
invariants PI, (PI would become an exact 2-forms. The 
adequate graphical language to account for all terms is 
the ribbon-like representation in which all plaquettes and 
Wilson contour are slided apart. The only troublesome 
contribution is the last one in Eq. (P| which explicitly 
depends upon the color orientation of the flux $((7) it¬ 
self. In the next section we analyze the arbitrariness of 
Sc and 7 -angles dependence of Eq. P|. 

D. Non-Abelian Bianchi Identities 

To complete the program outlined in sec. ITTT consider 
the surface independence of the non-Abelian Stokes theo¬ 
rem m- As one could expect the requirement of surface 
independence reduces to Eq. cni. On the other hand, 
the non-Abelian Stokes theorem PI applied formally to 
closed surface S'o gives 

^/(p)$(p) + fix = 2Trq{So), (30) 

pSSo xSSo 

where the integer q{So) is not vanishing in general and 
is discussed below. Since Eq. (PI is one of the central 
points of our work let us explicitly rederive it starting 
from Eqs. (0, p|. 

Consider Eq. o for some closed contour C 

U{CoC-^) = U{C) U{C-^) = U{C) U-\C) = 1, (31) 

part of which is shown on Figure^ There are two distinct 
surfaces Sc, Sq shown to the top and bottom of the con¬ 
tour with orientations induced by C. The non-Abelian 
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FIG. 4: Left: the non-Abelian Stokes theorem in application 
to the closed contour C which bounds two distinct surfaces 
Sc, S'c, SSc = SS'c = C. Arrows indicate the order of pla- 
quette vertices induced by the orientation of C. Right: the 
non-Abelian Bianchi identities for single lattice cube (see the 
text). 


have argued in sec. IIII Al Eq. (Tinil being applied to the 
infinitesimal cube is in fact the lattice implementation of 
the non-Abelian Bianchi identities and is illustrated on 
Figure 0 (right). It is clear that the integer q{So) is the 
magnetic charge discussed in sec. null Therefore, the 
non-Abelian Stokes theorem (EZI) which refers explicitly 
to the physically observable field strength allows to for¬ 
mulate the non-Abelian Bianchi identities on the lattice 
and to study their violation in gauge invariant terms. 


E. Discussions 


Stokes theorem (l?7ll applied for Sc and Sq leads to 

‘J’Sa = E HpMp)+J2 7x(^c)+2^ HSc) (32) 

pGSc xGSc x^C 

and analogous equation for d) 5 ^. The surface indepen¬ 
dence requires that $ 5 ^ = $ 5 ^ and therefore 

E ^ip)^(p) + E + E [ 7 : 1 :(-S'c) - 7a:(-S'c)] = 

P^Sq x^Sq x^C 

x^C 

= 2TT[k{S'c) - fc(5c)] = 2^ q{So). (33) 

Here Sq = Sc U and S'c is just the taken 
with reversed orientation due to which the terms 

changed sign in Eq. (TTHll . 
Consider the 7 -angles contribution in coming from 
points B,B' G C and let I^{ABC) denotes the Bargmann 
invariant for spinor wave functions at the points A, 
B, C. In particular, ^b{Sc) = A{A£B) and similarly for 
other 7 -angles. We note that one and the same unitary 
operator transforms A —> A', B ^ B', C ^ C . In other 
words the color directions of the fluxes at these points are 
rotated by one and the same rotation matrix. However, 
the Bargmann invariant being the area of the spherical 
triangle is unchanged when sphere is rotated. Therefore, 
the following identity holds 

A{ABC) - A{A'B'C') = 0. (34) 

It is clear that when Eq. m taken for each link of C is 
added to the l.h.s. of the total 7 -angles contribution 
becomes 

E[7:>:('5'c) - 7a:(-S'c)] = E 

x^C x^C 

where the orientation change of Sq in the inclusion 
:5o = 5c u S'c is crucial. For instance, VL[B) is given 
by A[A£'DC) and does not depend at all on contour C. 
We conclude therefore that Eq. (EOll is the consistency 
requirement for the non-Abelian Stokes theorem (I27II to 
be independent on the surface. But the point is that 
Eq. (IbOII is more than the consistency condition. As we 


This section is devoted to general notes concerning the 
Bianchi identities and the magnetic charge definition. We 
do not pretend on the exhaustive treatment, of course. 
However the following items seem to be worth mentioned: 

i) The SU(2) gauge invariance of the magnetic charge 
is evident from the fact that each term on the l.h.s. of 
CT is SU(2) gauge invariant by construction. The U(l) 
gauge invariance m of Eq. (IbOII is also obvious. One 
could argue that this Abelian symmetry is artificial and 
is only due to our intent to represent the non-Abelian 
flux direction in terms of the fictitious spinning particle. 
However, we do think that the U(l) invariance of IdOII 
might be relevant. Indeed, the interpretation of the Wil¬ 
son loop defining equation 03 in quantum mechanical 
language is natural and forces us to concentrate on the 
phase differences of wave functions (see, e.g., Eqs. flUll . 
(UHli, (IH), not on their concrete phases. Moreover, it 
allows to use the machinery related to the line bundle 
structure of quantum mechanics, mathematical founda¬ 
tions of geometrical phases and Bargmann invariants. In 
this respect the U(l) symmetry appears naturally and 
is inherent to our approach (it had been also discussed 
although in different context in Refs. 

ii) What was also crucial for our construction is the 
canonical orientation of elementary lattice plaquettes. 
We discussed this in details in sec. MB and concluded 
that in order to deduce the field strength from the in¬ 
finitesimal Wilson loops some canonical ordering must 
be introduced. It is true that in most cases the concrete 
ordering prescription does not matter since the usually 
considered quantities do not depend on it. For instance, 
the gauge action is insensitive to plaquette orientations, 
but this is certainly because the action is even in the field 
strength. As far as the magnetic fields are concerned 
their unambiguous definition is only possible with some 
canonical ordering prescription, overwise the components 
of could be determined only up to the sign even in 
the Abelian theory. However, it is clear that the ordering 
is not unique and although there are only few possibil¬ 
ities to choose from, the dependence of Eq. on the 
particular choice should be investigated separately. In 
this work we stuck with the conventional canonical or¬ 
dering described above, the ordering dependence will be 
investigated elsewhere. 
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Hi) As we have noted already it is natural to describe 
the non-Abelian Stokes theorem JTTIi in the ribbon-like 
graphical representation in which the theorem becomes 
essentially Abelian-like. In other words the non-Abelian 
nature of the theory is traded for the complicated geom¬ 
etry. Therefore the ribbon-like representation is actually 
not only the convenience. Once we could unambiguously 
assign each term in Eqs. f77|. ® to a particular geomet¬ 
rical object it is natural to ask whether these objects form 
a self-contained cell complex. For the non-Abelian Stokes 
theorem the answer is “no” because each Wilson contour 
requires the introduction of its own set of triangles (e.g., 
ABE on Figure 0 to which the 7 -angles are to be as¬ 
cribed. But the non-Abelian Bianchi identities do indeed 
allow the introduction of specihc cell complex in which 
every term on the l.h.s. of Eq. is unambiguously 
assigned to the particular 2-dimensional cell. Moreover, 
Eq. (HU could then be interpreted as usual coboundary 
operator acting on 2-cochains. Note that the above rea¬ 
soning resemble slightly the dual gravity-like representa¬ 
tion of SU(2) gluodynamics |3ISSEE 113 • We stress 

that this approach is not only the mathematical conve¬ 
nience. In fact it is the only way to analyze the structure 
of Eq. (Tinil at finite lattice spacing. In particular, it al¬ 
lows to show that the magnetic charge is closely related 
to the degenerate points 0 mentioned in sec. Ill Al f this is 
the topic of the next section). Here we note that the cell 
complex underlying Eq. (EOl) is described in Appendix 
the results of which are used in the next section. 

iv) It seems to be instructive to start from Eq. (EOJ, 
expand it in powers of the lattice spacing and get the 
Bianchi identities lU, lEJ in the continuum limit. How¬ 
ever, we failed to implement this program. As far as we 
can see the reason is two-fold. First, the original problem 
cu was posed quite differently from what could be ex¬ 
pected in the continuum. Indeed, our primary goal was 
to determine the magnetic charge and we intentionally re¬ 
fused to consider its gauge dependent color orientation. 
The manifestation of this could be seen by comparing 
Eqs. 1 ^, lI2Jl with (ll-iOII : while the former is in the adjoint 
representation and is vector in the color space the later is 
gauge invariant and is just one equation. Therefore it is a 
priori unclear how one could get 0,0 from 03) even in 
the limit of vanishing lattice spacing. On the other hand, 
Eq. H.Snil follows rigorously from 0 and we have no doubt 
that Eq. 03 indeed expresses the Bianchi identities on 
the lattice. Secondly, as we argue in the item below (see 
also the next section) the discussion of Eq. 13011 in the 
continuum limit is indispensable from the consideration 
of the degenerate points 0- 

v) Let us qualitatively consider what happens with the 
magnetic charge in the extreme weak coupling limit. 
The plaquette fluxes do not play any role since they are 
highly suppressed by the action. Therefore Eq. imui sim¬ 
plifies 

^ = 2nq{c), (36) 

xGSc 


where c is elementary lattice cube. Note that the mag¬ 
netic charge in not directly suppressed by the action and 
therefore there seems to be no reasons for it to die out in 
the continuum limit. Moreover, it is clear from ll^ that 
the non-zero q{c) is due to the particular distribution of 
the chromomagnetic field directions and is almost insen¬ 
sitive to the magnitude of the elementary fluxes. Indeed, 
each fix depends only on the flux directions and not on 
their magnitudes. In the next section we show that the 
non-zero r.h.s. of Eq. TO in the continuum limit indi¬ 
cates that at this point the chromomagnetic fields are 
degenerated and the particular determinant constructed 
from iff, Bf vanishes. 

IV. CHROMOMAGNETIC FIELDS 
DEGENERACY 

In this section we analyze the points of chromomagnetic 
fields degeneracy introduced in sec. ESI First we review 
the essential facts known in the continuum and then turn 
to the lattice dehnitions. 


A. Preliminaries 

In four dimensions the points of degeneracy of the chro¬ 
momagnetic fields are defined by 

detT = 0, (37) 


rjiab _ _ahc jpc _ ^ _ahc pc fQQ\ 

As we noted in sec. ITTO the physical signihcance of the 
points m crucially depends on the dimensionality. In¬ 
deed, in D=3 the operator coupled to the gauge poten¬ 
tials Af in the Bianchi identities 0 is 

rpab _ ^abc tdc _ ^ _ahc pc 

{3D) - £ (,39j 

and, in fact, is 3 x 9 matrix for which the determinant 
is undehned. We could at best consider the rank of the 
matrix TO and clearly 

rankr( 3 £i) < 3, (40) 

since is always the eigenvector with zero eigenvalue. 

We conclude therefore that in D=3 the very notion of 
chromomagnetic fields degeneracy is uncertain. 

In four dimensions the det T was calculated long 

ago [Illllil: 

detr oc det AT, (41) 


( 42 ) 
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ifoo = 2 detB, = = 

j^a _ rpa Tja _ 1 - jpa r\CL _ ^abc -rpb j^c 

~ 2 ^jk 1 Wk ~ ^ ^kij > 

where curly braces denote symmetrization B^^Qi^y = 
BiQk + BkQi- It is important that each element of K^i, 
is gauge invariant determinant constructed in terms of 
Ef and i?“. In particular, the off-diagonal elements are 
of the form det[ifiif/ci3/c], det[BiEkBk], where no sum¬ 
mation in k is implied and det[eie 2 e 3 ] is understood as 
the determinant of the column matrix constructed from 
color vectors 61 , 62 , 63 - Note that the off-diagonal ele¬ 
ments vanish identically in (anti)self-dual sectors [^ . 
However, our aim is not to analyze Eqs. CT-iHa in their 
generality. Rather we would like to show that the lattice 
Bianchi identities naturally lead to the same determi¬ 
nants ll^ . In particular, in the next section we show 
that the magnetic charge (HI is ultimately related to the 
zeros of these determinants and hence to the degenerate 
points m- 


B. detB = 0 on the Lattice 

In this section we consider first the three dimensional case 
which is much simpler geometrically. The results remain 
valid in four dimensions, but in D=4 there are important 
differences as well. 

It was noted in sec. rmn that the only reliable and 
rigorous way to analyze Eq. (isni) at finite lattice spac¬ 
ing is to consider the specially crafted cell complex for 
which Eq. pnji is the coboundary operation. The exis¬ 
tence and structure of this cell complex could be inferred 
by noting that the non-Abelian Stokes theorem and the 
Bianchi identities on the lattice are described most nat¬ 
urally in the ribbon-like graphical representation. Start¬ 
ing from this the cell complex could be completely con¬ 
structed (see Appendix). The advantage of this approach 
is that it is rather formal. Once we were able to assign 
a gauge invariant numbers (magnitudes of elementary 
fluxes and their relative orientations expressed in terms 
of Bargmann invariants) to each 2-dimensional cell, all 
we have to do is to consider the coboundary operator 
d : ^ C^, where is the fc-skeleton. For every lattice 

cube the action of d is equivalent to the Bianchi identities 
(EOl by construction and hence d assigns the correspond¬ 
ing magnetic charge to each lattice cube. However, the 
additional 0 -angles contribution implies that the geom¬ 
etry of the cell complex is not (hyper)cubical. In partic¬ 
ular, the 3-skeleton is larger than the union of lattice 
3-cubes. 

On the other hand, there is formally no difference be¬ 
tween different 3-cells of the complex. In particular, one 
can show that d : always assigns an integer 

number to every 3-cell. It is true that some of this “new” 
3-cells are trivial and the corresponding magnetic charge 
is always zero. However, there exist the non-trivial cases 
as well (see Appendix) one of which (and the only one in 




FIG. 5: Left: the flux directions in the plane (/i, i^) around 
point X in the weak coupling limit, eq lHHl . Right: the only 
non trivial 3-cell in D=3 (see the text). 


D=3) is illustrated on FigureEl (right). 

Consider some point x on the original D=3 lattice to¬ 
gether with 12 plaquettes and 8 cubes which share this 
point. Eq. applied to each cube forces us to take 
into account 8 triangles at cube’s corners (cf. Figure 0 
and to assign the corresponding Bargmann invariants Hi, 
i = I,..., 8 to each triangle. FigureH (right) shows the tri¬ 
angles around point x coming from different cubes. Note 
that all 8 triangles are properly oriented. By the same 
token one concludes that 6 squares, e.g. ABDE, are also 
valid 2 -cells of the cell complex and are equipped with the 
corresponding Bargmann invariants Ai(a;), i = I,. .., 6 . 
Then it is clear that the application of d : ^ to 

the set of 2-cells on Figure H assigns a well defined and 
gauge invariant integer number to the 3-cell shown on 
that figure: 


8 6 

27Tq{x) = + ^A^(x). (43) 

Formally it is just the same magnetic charge we have 
considered so far, but now it is ascribed to the site of 
the original lattice. We are confident that the magnetic 
charges in the lattice cubes correspond to the Bianchi 
identities violation. But what is violated in the lattice 
sites? 

To answer this question we expand m in powers of 
the lattice spacing. However, it is worth to mention that 
this expansion is not the usual one. In particular, it 
would be plainly wrong to look for O(a^) terms since the 
integer number on the l.h.s. of Eq. m does not depend 
at all on the lattice spacing. Therefore in the weak cou¬ 
pling expansion we should look for a-independent contri¬ 
butions or, better to say, to look for the conditions for 
a-independent terms to appear. 

In fact, all necessary relations were derived in Ref. [Tfi| . 
In particular, consider four plaquettes in the same plane 
which share the point x, Figure 0 (left). To the leading 
order the color directions of the fluxes at point x are 
given by 

X + fi + k : n-l-O(a^), 

X + fi — i) : ft — a D,^n + 0{a'^), 

X — jl + i) : fi — a Dfj^n + 0{a^), ^ ^ 

X — ft — u : n — a + D^)fi + 0{a^), 
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where we have denoted n = n(^^) (x) for brevity. We con¬ 
clude therefore that in the weak coupling limit the three 
points A,B,C (Figure 0 right) are distinguished: the 
flux directions assigned to them are in general indepen¬ 
dent and coincide with color direction of the particular 
component of The flux directions in all other ver¬ 
tices are obtainable by infinitesimal variation of the flux 
direction in one of the points A, B, C. 

Recall now that the Bargmann invariant assigned to 
each triangle and square is the oriented solid angle be¬ 
tween the corresponding flux directions. It follows then 
that the contribution of all squares is always of order 
0{a^) and is negligible. As far as the triangles are con¬ 
cerned they also give terms of order 0 {a^) unless the 
fluxes at points A, B, C become linearly dependent. In 
this case the corresponding Bargmann invariant could be 
±7r -I- 0{a) and the order 0{a) variation of the fluxes at 
various vertices is enough to change it by 2tt. It is clear 
that only in this degenerate case the non-zero l.h.s. of 
63 is at all possible. On the other hand, the flux direc¬ 
tions at the points A, C in the weak coupling limit are 
given by the corresponding chromomagnetic field compo¬ 
nents Bk- We conclude therefore that the non-vanishing 
magnetic charge 63 implies that the chromomagnetic 
fields are degenerate at this point 

q{x) 0 deti3(a:) = 0. (45) 

Note that the statement could not be reversed. For in¬ 
stance, in the case Bi = B 2 = B 3 both det B and q 
vanish. 

Eq. gni remains valid in four dimensions as well. The 
only distinction is that now we have 4 different mag¬ 
netic charges g^(a:) labeled by the direction fl dual to 
given 3-dimensional slice. In particular, the non-zero 
qfj.{x) implies that one of the determinants det[i3iB2S3], 
det[BiE 2 E 3 ], det [El i? 215 ^ 3 ], det[EiE 2 i? 3 ] vanishes. Note 
that these determinants are the diagonal entries of 
Eq. and therefore 

qfi{x) ^ 0 ^ = 0 (no sum over fi). (46) 

By symmetry considerations one expects that there 
should exist 3-cells for which the magnetic charge 
indicates the zeros of det[EiE 2 E 3 ], det[EiB 2 B 3 ], 
det[EiE 2 i? 3 ], det[BiB 2 B 3 ]. It turns out that these cells 
are (cc,/r, d^) (see Appendix). Indeed, the structure 
of cells is such that the argumentation leading to 
(im applies literally. Then the inspection of the flux di¬ 
rections assigned to vertices of shows that 

the non-zero magnetic charge of one of these 3-cells is the 
sufficient condition for the particular determinant above 
to vanish. 

As far as the off-diagonal elements of K^i, are con¬ 
cerned, they are highly sensitive to the topological 
properties of the gauge fields. Eor instance, K 12 = 
det[B 2 E 2 E 3 \ — det[EiEiE 3 ] vanishes in the (anti)self- 
dual sectors. It is possible to identify the 3-cells which are 
related to the off-diagonal entries of matrix. Indeed, 


consider the diamond-like 3-cells (see Appendix). In 
the weak coupling limit the flux directions assigned to 4 
plaquette corners become essentially the same and coin¬ 
cide with the corresponding component of . Then the 
flux orientations ascribed to 3 pairs of opposite vertices 
of are given by Efc, Efc, k = 1,2,3. Geometrically 
it is clear that for Ek = the 3-cells are highly 
degenerated and there is a good chance for the cobound¬ 
ary operator d : ^ ^( 2 ) give a non-zero magnetic 

charge. However, we are still lacking the rigorous argu¬ 
mentation here. One could only say (see also sec. ITO 
that the cells are indeed closely connected to the 
topological properties of the gauge background. The re¬ 
lation of the present approach to the gauge fields topol¬ 
ogy goes beyond the scope of the present publication and 
will be investigated elsewhere. 

To summarize, the non-Abelian Bianchi identities CT 
could be interpreted as the coboundary operator d : 

^ for the specific cell complex, the complicated 
geometry of which is the direct consequence of the non- 
Abelian nature of the theory. Moreover, the operator d 
considered in its generality necessitates the consideration 
of gauge invariant magnetic charges associated with var¬ 
ious 3-cells. While the non-vanishing magnetic charge in 
3-dimensional cube implies the violation of the Bianchi 
identities, in other 3-cells it is the sufficient condition 
for the particular determinant constructed from Ef,Bf 
to vanish. At finite lattice spacing these two types of 
magnetic charges are almost independent and should be 
considered as such especially since they are geometrically 
distinct: the former are ascribed to the lattice cubes, the 
later are assigned to the sites of the original lattice. How¬ 
ever, at vanishing lattice spacing the two types of mag¬ 
netic charges become closely interrelated (cf. Eq. (TIHll 'l: 
once the flux magnitude on the elementary plaquettes 
becomes negligible everywhere the non-Abelian Bianchi 
identities could only be violated at the degenerate points 

TO . 

V. NUMERICAL EXPERIMENTS 

It is true that the relevance of the above construction for 
the dynamics of the Yang-Mills fields is not evident from 
the preceding presentation. However, we specifically kept 
in mind from very beginning the possibility to apply our 
approach in real lattice experiments. In this section we 
describe the results of our numerical simulations. The 
problem to be considered is whether the violation of the 
Bianchi identities and the degeneracy of the chromomag¬ 
netic fields are physically significant. 

The general setup is as follows. We simulate the SU(2) 
lattice gluodynamics in three and four dimensions on the 
symmetric lattices with periodic boundary conditions. 
The action we adopt initially (see below) is the standard 
Wilson action. Until the sec. IvH the lattices we used 
are 16^ and 10^ with corresponding /3-ranges [5.0; 9.0] 
and [2.2; 2.8]. Note that these parameters are partially 
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unphysical. The purpose is to consider the behavior 
of the magnetic charges CT . (E51) . TO in various cir¬ 
cumstances, in particular, across the finite-volume phase 
transition. 

The simplest and instructive quantities to study are 
the densities p(/3), p(/3) of the magnetic charges (T^ . 
TO . The density p{j3) is defined irrespectively of the 
space-time dimensionality 

( 47 ) 

where summation is over all lattice 3-cubes and Nc is 
their total number. Evidently p(/3) measures the fraction 
of points at which the non-Abelian Bianchi identities are 
violated. The definition of p(/3) differs in D=3 and D=4. 
In three dimensions we have 

~P^^) = ^El«(")l- (48) 

where s is the lattice site, Ng is the total lattice vol¬ 
ume and q{s) was defined in sec. ITVBl In D=4 there are 
several types of the magnetic charges q and therefore the 
definition TO is ambiguous. We take the symmetric def¬ 
inition which looks similar to TO : s denotes the 3-cell 
which is not the lattice cube and Ng is the total number 
of these cells. Physically p{(3) is the fraction of the lat¬ 
tice volume occupied by zeros of various determinants, 
e.g. ll^ . (EHll . 

The dependence of p and p on /3 is shown on Figure El 
One can see that both densities are numerically similar in 
three and four dimensions and are almost /3-independent 
in accordance with general arguments of sec. Em In¬ 
deed, the /3-independence of p is certainly expected since 
there is no symmetry which could keep the sign of the 
determinants TO fixed. In particular, the perturbation 
theory gives the dominant contribution to the density 
p(/3). On the other hand the /3-independence of p fol¬ 
lows from the fact that the violation of Bianchi identities 
is closely related to the zeros of the above determinants. 
Therefore we come to the paradoxical conclusion that the 
perturbation theory also saturates the density p(/3). 

To resolve the problem we note that in the continuum 
limit the Bianchi identities are formulated for elemen¬ 
tary 3-volumes while the determinants are defined at any 
particular point. The corresponding construction on the 
lattice is essentially the same: the Bianchi identities and 
the magnetic charge q are ascribed to the elementary 3- 
cubes while the degeneracy points and the charge q are 
assigned to the lattice sites. It is important that these 
charges are geometrically distinct on the lattice: at arbi¬ 
trary small but non-zero spacing there is 0(a) distance 
between them and they are defined on different 3-cells. 
It turns out that on the lattice the magnetic charge at 3- 
cube and anti-charge at the neighboring site may coexist 
with almost no additional action penalty (cf. Eqs. (I3till . 
(133). Moreover, one can show that there could be no 
mechanism to prevent the creation of these ultraviolet 
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(UV) q-q pairs since it would violate the gauge invari¬ 
ance. Indeed, although the relative orientation of the 
fluxes is formally gauge invariant, any restriction of it 
will effectively squeeze the non-Abelian fluxes into one 
particular color direction. Then it would be hardly pos¬ 
sible to call the resulting theory non-Abelian |^. Note 
that the UV pairs above are irrelevant from the contin¬ 
uum viewpoint. Indeed, there is no trace whatsoever of 
the ultraviolet q-q pairs on the blocked lattice with lat¬ 
tice spacing N ■ a. At the same time the densities p(/3), 
p{P) account for all the charges q, q on equal footing and 
therefore are dominated by the UV fluctuations. 

We conclude therefore that the densities p(/3) and p(/3) 
are not the appropriate observables on the unblocked lat¬ 
tices. They are dominated by the ultraviolet noise which 
is only due to the mismatch in the domain of definition of 
the Bianchi identities and the degenerate points. It seems 
that the only way to make sense of the densities /?, p is 
to consider them on the blocked configurations for which 
the ultraviolet noise is gradually removed. However, our 
approach to the problem is different and is described be¬ 
low. 


A. Modiflcation of the Action 

As follows from the above presentation, the dynamics of 
q and q magnetic charges is highly UV sensitive and the 
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dominant configurations are small (at the scale of UV 
cutoff) q-q pairs. It seems that this observation forbids 
the discussion of the significance of the Bianchi identi¬ 
ties violation and the points of degeneracy since it is im¬ 
possible to separate the UV noise from physically rele¬ 
vant excitations. Essentially the same problem exists in 
usual field theories, where the vacuum condensates are 
commonly used to parametrize the non-perturbative ef¬ 
fects. The well known example is the gluon condensate 
(as(E“^)^) which is perturbatively divergent but its non- 
perturbative part is non-vanishing and is known to be of 
major phenomenological importance. The subtraction of 
the perturbative tail of various condensates is challeng¬ 
ing and the usual approach is to subtract it order by 
order in the coupling constant. However, we don’t see 
any tractable way to do this in our case. 

On the other hand, it is possible to reformulate slightly 
the original problem. Instead of trying to isolate the 
effects due to the UV q-q pairs we could equally ask 
what happens when the magnetic charges are partially 
removed from the vacuum. Indeed, the definition of q 
and q charges is local and gauge invariant. Therefore, 
nothing prevents us from modifying the Wilson action to 
include the additional terms which could influence the 
dynamics of q, q charges. Since it is hardly possible to 
invent the additional well defined terms which are sensi¬ 
tive to the UV dynamics only, we will study the following 
simplest modification 

p C S 

where the first term is the standard Wilson action and 
c denotes the elementary lattice cubes. The last term 
in Eq. has different interpretation in three and four 
dimensions. In D=3 s denotes the lattice sites and q{s) 
is given by Eq. 63. In four dimensions the last term 
a priori depends on the concrete definition of the mag¬ 
netic charges q. As in the previous section we take the 
symmetric definition: s denotes the 3-cells which are not 
the lattice cubes and g(s) is the corresponding magnetic 
charge. It turns out that our results are almost insensi¬ 
tive to the particular choice of the last term in Eq. (03, 
see sec. IV (II 

The modified action is local and SU(2) gauge invari¬ 
ant. Indeed, from the defining equations CT . 03 one 
can see that 03 intertwines the links which are at most 
two lattice spacings apart, while the gauge invariance 
follows by construction. Then the universality suggests 
that the continuum limit of the model defined by 03 
should be the same as one for the model with the con¬ 
ventional Wilson action (see also sec. I VII for discussions). 
On the other hand, the additional coupling constants 7 , 
7 allow to study the effects which are due to the Bianchi 
identities violations and the degeneracy points. The par¬ 
ticular limit 7 —!■ 00 is of special interest since it cor¬ 
responds to the theory with nowhere violated Bianchi 
identities. As far as the 7 coupling is concerned we are 
not so confident that the limit 7 ^ 00 corresponds to a 


sensible theory. For instance, in D=3 the nowhere van¬ 
ishing det B = det[BiB 2 B^] implies that it is of the same 
sign everywhere, which contradicts the perturbative ex¬ 
pectations [U and probably violates CP symmetry. At 
the same time the point 7 = 7 = 0 is certainly equivalent 
to the conventional lattice gluodynamics. 

In the next two sections we study the model 03 along 
the lines 7 = 0 and 7 = 0 in the ( 7 , 7 ) parameter space at 
fixed value of the gauge coupling /3. The simulations were 
performed on 20^ and 12^ lattices at /3 = 6.0 and /3 = 2.4 
correspondingly. Note that this choice of parameters is 
based on the experience with pure YM theory, in which 
these /3 values and volumes correspond to the physical 
scaling regime . While the point 7 = 7 = 0 was 

simulated with standard overrelaxed heatbath updating, 
away from it we implemented the Metropolis algorithm 
which is the only one available at non-zero 7 , 7 . The 
procedure turns out to be very time consuming especially 
in D=4. Indeed, the one link update step requires to take 
into account the magnetic charges q, q in all neighboring 
cells the number of which is much larger in D=4 (see 
Appendix). Because of this we were unable to thoroughly 
scan the ample range of 7 -couplings, only the following 
points were considered in details 

( 7 . 7 ) 3 D = {(0,0); (4,0), (7,0), (9,0); (0,4)}, 

( 7 . 7 ) 4 D = 1(0,0); (4,0), (6,0), (8,0); (0,4)}. 

In particular, the complexity of the algorithm precludes 
us from studying the phase diagram of the model 03 
(see below) and investigate the finite volume effects. Be¬ 
low it is silently assumed that the chosen volumes are 
large enough even at non-zero 7 , 7 couplings. At each 
7 -point we generated about one hundred statistically in¬ 
dependent gauge samples separated by ~ 10^ Monte 
Carlo sweeps. The observables of primary importance 
are the planar Wilson loops from which we extracted the 
heavy quark potential (see, e.g. Ref. for details) 
and the correlator of the Polyakov lines {P{0)P{R)), 
P{x) = 1/2 Tr}}^ t7o(x -I- t). To improve the statistics 
the standard spatial smearing and hypercubic block¬ 
ing for temporal links were used. In D=4 we also 
monitored the topological charge Q, the topological sus¬ 
ceptibility X = (Q ^ )/V defined by means of the overlap 
Dirac operator (see, e.g. Ref for details and 
further references). 


B. 7 = 0 Line 

Here we study the effect of the gradual removal of the 
points in which the Bianchi identities are violated. Let 
us consider first the behavior of the densities 03 , 03 
with raising 7 coupling. It turns out that p stays al¬ 
most constant (Figure |3 upper panel) in both three and 
four dimensions monotonically varying from 0.477(1) to 
0.435(1) in D=3 and from 0.2831(1) to 0.2441(2) in D=4 
in the entire 7 -range considered. On the other hand. 
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the density p falls down exponentially with 7 and be¬ 
comes of order 0(10“^) in D=3 (0(10“^) in D=4). We 
note in passing that the mean plaquette (l/ 2 Trt/p) is 
also almost insensitive to the 7 coupling (Figure [ 7 | bot¬ 
tom) rising in D=3 from 0.8248(1) to 0.8263(3) when 7 
is changed in the entire range (the corresponding change 
in D=4 is from 0.6301(2) to 0.6548(2)). 

A few comments are now in order. First, the con¬ 
stancy of p and the simultaneous falloff of p by orders 
of magnitude implies that the above picture of dominat¬ 
ing ultraviolet q-q pairs is greatly oversimplified. It seems 
that the UV fluctuations are indeed dominating, but their 
structure is much more involved. In particular, it has lit¬ 
tle to do with the model of tightly bounded q-q dipoles, 
rather it is some complicated mixture of various charge- 
anticharge configurations which probably don’t form the 
dipole-like pairs at all. 

Turn now to the behavior of the heavy quark potential 
and the Polyakov lines correlation function with rising 7 
coupling. In three dimensions (Figure |HJ) both the Wil¬ 
son loops and the correlator (P(O)P(i?)) show almost no 
sign of 7 coupling dependence, in particular, the asymp¬ 
totic string tension at large 7 is equal to its value in the 
pure Yang-Mills theory. However, the situation changes 
drastically in D=4. One can see from Figure |21 that the 
correlation function (P(O)P(i?)) tends to non-zero posi¬ 
tive value at large separations when 7 coupling becomes 
of order few units 

lim (P(0)P(i?)).v>i = const >0. (51) 

R—>^oo ~ 

The heavy quark potential extracted from Wilson loops 
is shown on Figure HOI (upper panel) and for 7 ^ 1 is 
indeed flattening at large distances 

lim Vy>i{R) = const. (52) 

Note that it is hardly possible to conclude firmly from 
Figure im along that the asymptotic string tension is in¬ 
deed vanishing; however, Eq. 15111 and Figure IHI are in¬ 
compatible with its non-zero value. 

The other measured observables do not show strong 
dependence on 7 coupling. In particular, the topological 
charge Q stays at zero in average value albeit with slightly 
narrower distribution. As is clear from the bottom panel 
of Figure ^ the topological susceptibility x = {Q'^)/V 
diminishes at 7 « 1 by approximately 25% and the esti¬ 
mation of its limiting value is 

lim x^^'^( 7 ) = 163(8) MeV, (53) 

7—^00 

which should be compared with = 

212(3)MeV in pure YM theory, where the physical units 
are fixed by the string tension = 440 MeV. 

The discussion of the results presented above is post¬ 
poned until sec. El Here we only note that the dynam¬ 
ics of YM fields in D=3 seems to be almost insensitive to 
whether or not the Bianchi identities are violated. In par¬ 
ticular, the complete suppression of the magnetic charges 


which indicates the violation of the Bianchi identities has 
almost no consequences for the correlators we considered. 
However, the four dimensional case appears to be quite 
different. Our results indicate that the suppression of the 
Bianchi identities violation is likely to destroy confine¬ 
ment while other measured characteristics of the theory 
remain essentially unchanged. 


C. Suppressing the Degenerate Points 

Consider the response of the theory on the suppression 
of the degenerate points. The qualitative difference in 
the behavior of the system along the lines 7 = 0 and 
7 = 0 could be seen already on the simplest observables 
like p, p. We have checked that the falloff of the degen¬ 
erate points fraction piff) is indeed exponential with 7 
in both three and four dimensions; the relevant num¬ 
bers are P 3 _d( 0 ) = 0.477(1), P 3 _d( 4 ) = 0.0098(2) and 
P 4 ^( 0 ) = 0.2831(1), P 4 d( 4 ) = 0.045(3). However, the 
fraction of points at which the Bianchi identities are vi¬ 
olated also notably diminishes with 7 . The falloff of 
p( 7 ) in D=3 is not so pronounced (p(0) = 0.1758(4), 
/o(4) = 0.1010(4)) and starting from 7 « 1 it is numeri¬ 
cally larger than p. It is surprising, however, that in D=4 
the inequality p < p holds for all 7 values considered and 
in fact the fraction of points at which the Bianchi identi¬ 
ties are violated is diminished by the order of magnitude 
(0.2059(2) at 7 = 0 versus 0.024(3) at 7 = 4). As far as 
the mean plaquette energy is concerned its behavior is 
similar to that on the 7 = 0 line. In particular, in D=3 
it essentially stays constant while in four dimensions it 
changes from 0.6301(2) to 0.655(1). 

As we noted already the suppression of the degenerate 
points ll^ . (EH) might not be physically meaningful. For 
instance, in three dimensions the orientation of the triple 
(Hi, i? 2 , H 3 ), although being gauge invariant, is not fixed 
by any symmetry or physical principle. The attempt to 
fix the sign of det B everywhere probably will lead to 
physically unacceptable results. Indeed, the closer in¬ 
spection of the Polyakov lines correlator reveals that it is 
an oscillating function of the distance. Hence the lattice 
reflection positivity is lost and the theory seems to be 
pathological at non-zero 7 . 

In four dimensions the suppression of the degenerate 
points leads to qualitatively the same results which, how¬ 
ever, are much more pronounced. For instance, the Wil¬ 
son loops (yV{R,T)) measured at 7 ^ 0 are notably os¬ 
cillating at fixed R and varying T iFigureim top panel). 
However, unlike the three dimensional case we can easily 
pin-point the origin of the reflection positivity violation. 
Indeed, it is well known that in the fixed topological sec¬ 
tor the theory certainly violates CP and it is natural 
then to ask what is the typical topological charge of the 
configurations at non-zero 7 . 

The bottom panel of the Figure ^2 shows the Monte 
Carlo history of the topological charge on 8 ^ lattice at 
f3 = 2.30, 7 = 0, 7 = 0.1, 0.5 when the starting configura- 
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tion was thermalized at 7 = 7 = 0 (note that we changed 
the lattice geometry and the /3 coupling for reasons to 
be explained shortly). In view of the observed reflection 
positivity violation at 7 7 ^ 0 it is not surprising that Q in¬ 
deed stays away from zero in average. What is surprising, 
however, is that the average topological charge {Q) turns 
out to be always positive and extremely large for 7 > 0 . 
In particular, for 0 < 7 XC 1 the mean topological charge 
is shifted only slightly from zero being of order few units. 
However, once the 7 coupling becomes comparable with 
unity Q flows away from zero during Monte Carlo updat¬ 
ing towards extremely large positive values with almost 
constant and very high rate. In fact, it quickly becomes 
too large to be technically accessible for us and this was 
essentially the reason to consider so small lattices here. 
The volume dependence of (Q) could be inferred by not¬ 
ing that the last term in ( 1 ^ responsible for the rapid 
increase of the topological charge is the bulk quantity. 
Therefore {Q) seems to be proportional to the volume 
at fixed 7 although we had not thoroughly investigated 
this dependence numerically. We have checked that the 
behavior of Q at non-zero 7 is always similar to that on 
Figure m irrespectively to the concrete meaning of the 
last term in the action It does not matter which 

particular type of magnetic charges q is suppressed by 
the 7 coupling, we always see the violation of the reflec¬ 
tion positivity which is due to the rapid increase of the 
global topological charge. This problem is discussed in 


next section. 


VI. DISCUSSIONS 

The interpretation of the results we achieved so far may 
not be simple and straightforward. Here we discuss a few 
particular points which are essential for our work. 

First of all, we do see that the physical significance of 
the Bianchi identities is quite different in D=3 and D=4. 
The three dimensional theory turns out to be insensi¬ 
tive to the suppression of the Bianchi identities violation. 
Even the complete removal of q charges from the vacuum 
does not change the theory in any notable way. The four 
dimensional theory seems to be different in this respect. 

The suppression of the Bianchi identities violation is 
likely to destroy confinement liberating color charges in 
the fundamental representation. It is tempting to con¬ 
clude then that the confinement phenomenon is due to 
the field configurations for which the r.h.s. of Eq. m is 
non-vanishing. This conclusion looks natural for the fol¬ 
lowing reasons. First, it matches the known confinement 
mechanism in the simple Abelian models. Secondly, it 
could explain why in the continuum considerations con¬ 
finement is missing since usually the Bianchi identities 
with vanishing r.h.s. are taken for granted. Third, it 
qualitatively matches the phenomenological lattice ob¬ 
servations that the geometrically thin line- or string-like 
objects (Abelian monopoles, P-vortices) might be rele¬ 
vant for confinement (see, e.g Refs. for review and 

further references). And finally it does not look hopeless 
from the field-theoretical point of view since, as we ar¬ 
gued above, at vanishing lattice spacing the mechanism 
of the Bianchi identities violation has little to do with 
singular fields, rather it is related in some complicated 
way to the points of chromomagnetic fields degeneracy. 

However, the striking difference between three and four 
dimensional theories with respect to the suppression of 
the Bianchi identities violation shows that this conclusion 
is probably misleading. If the confinement phenomenon 
is indeed due to the Bianchi identities violation then it 
should disappear also in D=3 at large 7 coupling. But 
this does not happen and hence we come to the unnatural 
conclusion that the confinement mechanism has little in 
common in D=3 and D=4. 

However, we could take a different point of view. 
Namely, there is indeed a great physical difference be¬ 
tween the Bianchi identities in three and four dimen¬ 
sions. As we discussed in sec. EH Eq. o constitutes the 
algebraic restriction on the gauge potentials for a given 
distribution of the chromomagnetic fields. Away from 
the degeneracy points [s^ the gauge potentials could be 
completely reconstructed just from the Bianchi identi¬ 
ties along. In this respect the violation of the Bianchi 
identities could be seen as the source of the gauge poten¬ 
tials ambiguities and the suppression of non-zero r.h.s. 
of Eq. (Q effectively restricts the gauge inequivalent A“ 
which are to be taken into account in the functional inte- 
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gral. It is crucial that in three dimensions the analogous 
argumentation fails and in fact Eq. does not restrict 
the gauge potentials in any notable way irrespectively 
whether or not it is violated. 

The natural and probably the only available quantity 
which is sensitive to the gauge potentials ambiguities is 
the (^min) condensate. Therefore the following qualita¬ 
tive scenario emerges. It is the non-perturbative 
condensate which seems to be relevant for confinement. 
In four dimensions the Bianchi identities are the tool 
which allows to restrict the (^^m) condensate. More¬ 
over, the suppression of non-zero r.h.s. of Eq. o makes 
to vanish. Clearly the same approach does not 
work in three dimensions because the Bianchi identities 
do not constraint in D=3. Note that this is only 
the qualitative picture. In particular, the dependence of 
on the q charges density could be very compli¬ 
cated especially because of the dominating perturbative 
contributions. 

To reiterate the point we note that the justification 
to consider the modified action is that the second 
term in ( 1 ^ is local and preserves all the symmetries of 
the original action. Then the universality suggests that 
the continuum limit of the model (B51) should be 7 cou¬ 
pling independent (we take 7 = 0 for definiteness). At 
the same time our results indicate that this is probably 
not the case. If we would accept the Bianchi identities 
violation as the primary reason for confinement then we 
would be faced with serious universality problems. How¬ 
ever, the above scenario based on the condensate 

seems to avoid (at least formally) this issue. 

In fact, the dependence of (A^j^j) condensate on the 7 
coupling could be measured directly. Namely, we could 
measure the quantity (A^) in the Landau gauge and its 
drop with rising 7 coupling gives an estimate for the be¬ 
havior of the (A^j„) condensate when the violation of 
the Bianchi identities is gradually removed. Moreover, 
this could be compared with the results of Ref. @] where 
the same Landau gauge (A^) albeit with different nor¬ 
malization was measured across the finite-temperature 
deconfinement phase transition. Note that the quantity 
(A 2 ) in the Landau and Coulomb gauges was already 
introduced in Refs, The details of our measure¬ 

ments are as follows. The gauge potentials are defined in 
terms of the link matrices U^{x) 

K = Tr|l[C/4a:)-[/t(x)], (54) 

where a is the lattice spacing. The Landau gauge was 
fixed by minimizing with overrelaxation 

algorithm until the magnitude of 9^A“ becomes every¬ 
where less than I0“®. The results are presented on Fig¬ 
ure ESI One can see that in three dimensions (A^) is in¬ 
deed almost insensitive to the 7 coupling confirming the 
qualitative scenario outlined above. On the other hand, 
in four dimensions (A^) drops down with increasing 7 
by essentially the same amount which was reported in 
Ref. [2- Note that the relative drop of (A^) is expected 



FIG. 12: The quantity (A^) in the Landau gauge at 7 = 0 as 
function of 7 in three and four dimensions. 

to be small nil Indeed, on general grounds we have 



and clearly the Landau gauge (A^) is dominated by the 
perturbative tail at weak coupling. However, the drop in 
(A^) across the phase transition is believed to be entirely 
due to the non-perturbative condensate (A^^^). 

The next comment concerns the behavior of the topo¬ 
logical charge with respect to the degenerate points sup¬ 
pression. It is true that the violation of the reflection 
positivity with rising 7 coupling is to be expected on the 
general grounds. Moreover, it is also expected that in 
D=4 the non-zero in average global topological charge 
is the origin of the reflection positivity violation. How¬ 
ever, the following questions remain: why the topological 
charge in always positive and rises so rapidly with respect 
to the Monte Carlo updating? What is the relation be¬ 
tween Q and the magnetic charges g, ql 

Evidently, the non-zero in average value of the topo¬ 
logical charge requires it to be mostly either positive or 
negative. At the same time our experience with various 
possible definitions of the last term in the action 
shows that the positivity of Q at 7 0 is seemingly 

built in our approach from very beginning. As far as 
we can see the only place which distinguishes between 
<5 ^ 0 is the canonical orientation of the elementary pla- 
quettes which we accepted (see sec. irrmii . Indeed, the 
construction of the magnetic charges g, q depends on the 
particular canonical orientation which is not uniquely de¬ 
fined. Although there are only few possibilities to choose 
from, different choices could discriminate the sign of the 
topological charge. Indeed, in the fermionic language the 
sign of Q distinguishes the left and right chiralities (ori¬ 
entations) analogously to the canonical orientation which 
discriminates the left and right coordinate systems. Thus 
we expect that the sign of Q at non-zero 7 will change 
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with inequivalent choice of the canonical orientation. As 
far as the rapid growth of Q is concerned it seems that 
the only possible explanation is that the non-zero 7 cou¬ 
pling lifts the degeneracy of different topological sectors. 
The situation is reminiscent to the quantum mechanical 
problem of the periodic potential on which the constant 
7 -dependent electric field is superimposed. 

It would be instructive to have the explicit expression 
for Q in terms of the q, q magnetic charges distribution. 
However, the usual approaches available in the literature 
(see, e.g. Refs, [^l^) seemingly lead to the erroneous 
results. For instance, the treatment of Ref. applies 
almost literally in our case. The outcome is that the topo¬ 
logical charge is given by the linear combination of the 
q, q magnetic charges and hence vanishes when g, q are 
highly suppressed. This seems to contradict the observed 
rapid growth of Q with rising 7,7 couplings when all the 
q, q charges are suppressed on equal footing. Therefore 
either the results of Ref. [t^ should be modified in our 
case or we should look for different definition of the topo¬ 
logical charge. The definition of the topological charge, 
which closely follows the approach of the present work, 
could be given along the lines of Refs. Izlizi (see also 
Ref. [l^ for an excellent introduction). Instead of study¬ 
ing the evolution of single spinor along the closed con¬ 
tours we could consider the corresponding evolution of 
the degenerate two-level system for which the resulting 
non-Abelian geometrical phase describes the YM instan- 
ton. This approach is under investigation and will be 
published elsewhere. 


VII. CONCLUSIONS 

In this paper we considered the non-Abelian Bianchi 
identities in SU(2) pure Yang-Mills theory focusing on 
the physical significance of the chromomagnetic fields de¬ 
generacy points and the possibility of Bianchi identities 
violation. These questions necessitate the regularization 
and we specifically kept in mind the lattice formulation. 
It had been known for a long time that the Bianchi identi¬ 
ties in general are the requirement that the gauge holon- 
omy for any null-homotopic path equals to unity. The 
main achievement of this paper is the reformulation of 
the above requirement in terms of the physical elemen¬ 
tary fluxes (field strength). Our approach is based on 
the non-Abelian Stokes theorem appeared recently and 
allows to give an explicit gauge invariant expression for 
the Bianchi identities on the lattice. Simultaneously it al¬ 
lows to formulate the notion of the non-Abelian Bianchi 
identities violation in gauge invariant and local form. 

As a further development of our approach we showed 
that the study of the lattice Bianchi identities natu¬ 
rally leads to the consideration of the chromomagnetic 
fields degeneracy points at which a particular determi¬ 
nants constructed from Ef, Bf vanish. It turns out that 
the violation of the Bianchi identities and the degenerate 
points are closely related to each other. In particular, in 


the weak coupling regime the Bianchi identities violation 
is not related generically to the singular fields, rather it 
is due to the existence of the degenerate points. 

As is clear from the above presentation the main ad¬ 
vantage of our approach is that the non-Abelian nature of 
the theory had been traded for the complicated geometry 
which, however, allows the pure geometrical Abelian-like 
treatment. Then both the Bianchi identities violation 
and the degeneracy points formally appear as usual mag¬ 
netic charges. However, we stress that the term “mag¬ 
netic charge” and, in fact, the entire Abelian analogy is 
only formal. In particular, the physical interpretation 
oi q, q charges is completely different; there is no mag¬ 
netic charge conservation whatsoever on the original (hy- 
per)cubical lattice. Nevertheless, the Abelian-like repre¬ 
sentation is invaluable for the analysis presented above. 

The locality and gauge invariance of the definition of 
the Bianchi identities and the chromomagnetic fields de¬ 
generacy points permits us to modify the original gauge 
action and to study the effects of gradual removal of these 
objects from the vacuum. It turns out that in the four 
dimensional case the suppression of the Bianchi identi¬ 
ties violation seems to be relevant for confinement: the 
heavy quark potential extracted from Wilson loops flat¬ 
tens at large distances and the correlator of the Polyakov 
lines tends to non-zero constant at large separations. At 
least this is the case on the lattices we have studied. At 
the same time, other correlation functions which we mea¬ 
sured had not been changed considerably. The situation 
in D=3 turns out to be just opposite. Namely, the the¬ 
ory is almost insensitive to the suppression of the Bianchi 
identities violation. However, in D=4 the complexity of 
the numerical simulations precluded us from studying the 
relevant issues like the phase diagram of the modified 
model, the volume dependence of our results, etc. We 
hope to address these questions elsewhere. 

As far as the degenerate points are concerned any at¬ 
tempt to remove them from the vacuum results in the 
reflection positivity violation. Moreover, in D=4 this vi¬ 
olation is due to the extremely large positive global topo¬ 
logical charge which grows rapidly during Monte Carlo 
updating. This observation could be relevant for study¬ 
ing the gluodynamics in the topologically non-trivial sec¬ 
tors. 

Confronting the results obtained in D=3,4 we argued 
that it is probably misleading to consider the violation 
of the Bianchi identities as the primary cause of con¬ 
finement. Instead the correct picture would be to inter¬ 
pret the Bianchi identities as an algebraic constraint on 
the gauge potentials and to relate the confinement phe¬ 
nomenon to the existence of the non-perturbative (A^j„) 
condensate. This scenario seems to be in agreement with 
universality expectations, works the same in both three 
and four dimensions and does not contradict our findings. 


19 


Acknowledgments 

The authors are grateful to prof. V.I. Zakharov, 
E.T. Akhmedov and to the members of ITEP lat¬ 
tice group for stimulating discussions. This work 
was partially supported by grants RFBR-03-02-16941, 
RFBR-05-02-16306a, RFBR-05-02-17642 and RFBR- 
0402-16079. The work of F.V.G. was supported by IN- 
TAS YS grant 04-83-3943. 


APPENDIX 


Here we describe the cell complex underlying the Bianchi 
identities CT . We start from single plaquette and note 
that the application of non-Abelian Stokes theorem 
assigns spinor wave function (z | to each plaquette cor¬ 
ner. This could be represented by 4 points belonging to 
this plaquette and shifted from the corners towards the 
plaquette center. The totality of this points constitutes 
the 0 -skeleton C° of the cell complex and it is convenient 
to parametrize s G by the point x of the original lat¬ 
tice and by two shifts with corresponding shift directions 
(see Figure El left): 


s = s(x, n, V,dv) = s{x, V, dy, n,d^), 
fx ^ ly ^ di/ — zbl. 


(55) 


In total there are 2D{D — 1) ■ V sites, where V is the 
lattice volume. 

Turn now to the 1-skeleton which consists of two 
types of links. The first group contains the original links 


Si = s(x, /i, d^, ly, du) ^ Sf = s{x + fi, fi, -d^, V, d„), (56) 

which carry the matrix element ( z(si) | I7^(x) | z(s/)). 
Links from the second group 


Si = s{x,n,d^,iy, d„) 


^ ^ A 


Sf = s{x,\,dx,v,dy) 


(57) 

are ascribed with the matrix element ( z{si) \ z{sf)). 

As far as the 2 -skeleton is concerned its structure is 
different in three and four dimensions. As a consequence 
C^, k > 2 , also differ considerably and are described 
separately below. 


constructed from links m- Namely, (x,/i, d^) has 
the following four points in its boundary 

y^^\x,n,df^) : s{x, n,df_t, v,±d^), s(a:, d^. A, ±dA), 

II ^ V ^ \ . (58) 

The coboundary operator d : ^ assigns the 

Bargmann invariant (EHli to the particular square. The 
argumentation of sec. mm allows to show that 

1) = , (59) 

where we have denoted the values assigned to each square 
by the same symbol hoping that this will not lead 

to confusion. Third type of 2-cells contains various tri¬ 
angles constructed from links (15711 : there are three 
points in the boundary of 

: six, n,d^, i/,d^), (60) 

six, V, du. A, dx), six, A, dx, d^), 
fi ^ v ^ X . 

The operator d : ^ assigns the corresponding 

Bargmann invariant to the triangle. Note that the last 
group of 2-cells is formed by mixture of links (ISBll . itFzl) 
and need not be considered, in fact: by Eq. the phase 
associated with them is always zero. It is important that 
the value assigned by d to every 2-cell is always taken 
modulo 27r and is rather similar to Opiaq = [dd] 27 r in the 
language of compact U(l) gauge model. In other words 
it is silently assumed that only gauge invariant quantities 
are ascribed to every 2-cell. 

As far as the 3-skeleton is concerned it contains 
essentially two types of 3-cells. First, the original lat¬ 
tice cubes which look as on Figure 01 (right); each cube 
contains 6 plaquettes and 8 triangles at its corners. The 
coboundary operator d : C'^ —>■ considered for any 
particular cube is identical to Eq. by construction. 
The 3-cells of the second group are constructed entirely 
from triangles and squares above and are illus¬ 
trated on Figure El (right). The physical meaning of the 
corresponding magnetic charge is analyzed in sec. nm 
Thus the consideration of three dimensional case is 
completed. Note that geometrically there is one more 
type of 3-cells which, however, need not be taken into 
account. These 3-cells are formed by two squares (EHl) 
and four links connecting them. It follows from 
and that d : ^ always gives zero on these 

cells. 


D=3 

Here contains three types of 2-cells. First, there are 
original plaquettes p the boundary of which consists of 
the links (EH). Moreover, the standard coboundary op¬ 
erator d : ^ acts in accordance with Eq. (E51) and 

assigns the flux magnitude 4>(p) to the plaquette. Second 
group of 2 -cells contains various squares ix, p,dp) 


D=4 

In four dimensions the consideration of the cell complex 
underlying the lattice Bianchi identities (EOl) becomes 
cumbersome. In particular, we do not give the full list 
of cells forming C^, k = 2,3, 4, only cells relevant to the 
considerations in sec. II V BllVl are presented. 
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FIG. 13: 3-cell in D=4 (see the text). 

First we note that the D=3 construction applies di¬ 
rectly in D=4. In particular, 2-skeleton includes the pla- 
quettes, squares JHHl and triangles (HU, trivially general¬ 
ized to four dimensions. In the 3-skeleton we identify 
then the usual 3-cubes and 3-cells shown on Figure 0 
(right). 

However, it is clear that in D=4 the C^, are not 
exhausted by the above 2- and 3-cells. In particular, 
the 2-skeleton contains now an additional set of triangles 
with vertices 

s{x,fi,d^,u,d^), (61) 

s{x,fj,,df_,,X,d\), s{x, fj.,dfj., p,dp), 

p, ^ V ^ \ ^ p, 

and squares , the vertices of which are 

: s{x,p,dfj,,u,d^), s{x,iy,d^,X,dx), (62) 

s{x,X,dx,p,dp), s{x,p,dp,p,dp), 
p ^ V ^ X ^ p. 


All these 2-cells are constructed from links iZJ and 
therefore are ascribed with the appropriate Bargmann 
invariants. 

In the 3-skeleton the new diamond-like cells ^ con¬ 
sisting of 6 vertices and 8 triangles appear. In turn, these 
3-cells could be subdivided into two groups. 

{x, p,dp): the 3-cells in this group are constructed 
from 8 triangles lICTl) and are similar to those considered 
in D=3, Figure H (right). In particular, one can show 
that (cf. Eq. 1^111) 

^^^\x,p,l) = ^^^\x + p, p,-1), (63) 

see the note following Eq. HU . The physical interpreta¬ 
tion of the corresponding magnetic charge is discussed in 
sec. lIVBl 

(a:, dp, dy, d\, dp): these 3-cells are built from both 
types of triangles (Rinil . lICTll . The corresponding vertices 
are constructed by fixing a particular combination of shift 
directions dp-, there are 6 distinct planes passing through 
given lattice site in D=4 and s(x, p,dp,v,du), p ^ v \s 
one of the six vertices of dp, d^, d\, dp) cell. The 

total number of these 3-cells per lattice site is 2"* = 16. 
Note the specific pattern of the flux directions assigned 
to the vertices of {x, dp,dv,d\, dp), which is radically 
different from what we have encountered so far. In the 
weak coupling limit the opposite vertices are ascribed 
with the same components of chromoelectric and chro- 
momagnetic fields (see the Figure ICTl. 
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